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We report on recent numerical computations of the Landau gauge gluon and ghost propagators as well as of the 
ghost-gluon-vertex function in pure SU (3) Yang-Mills theory and in full QCD on the lattice. Special emphasis is 
paid to the low momentum region. In particular, we present new data for the gluon propagator at momenta below 
300 MeV. We also discuss different systematic effects as there are finite-size, lattice discretization and Gribov 
copy but also unquenching effects. A MOM-scheme running coupling CXj(^^) based on the ghost-gluon vertex 
is calculated and found to decrease for momenta below 550 MeV, even though the renormalization constant of 
the vertex deviates only weakly from being constant. 
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I. INTRODUCTION 



With our contribution to this workshop we give an overview 
on recent lattice computations of the Landau gauge gluon and 
ghost propagators in quenched (Nf — 0) and in full QCD 
(Nf = 2). The full analysis and additional results of related 
observables can be found in the Ph.D. thesis of one of us |[ll]. 

For QCD being the theory of the strong interaction a co- 
herent description of all hadronic features directly based on 
the dynamics of confined quarks and gluons, given in terms of 
all propagators and vertex functions of QCD, should be avail- 
able. As reported by R. Alkofer at this conference, they may 
serve as an input from first principles for the Bethe-Salpeter 
and Faddeev equations and this opens a way to a model inde- 
pendent phenomenology of nonperturbative phenomena. 

Lattice computations of gauge-variant Green functions 
have attracted more and more interest in recent years, because 
the results can be directly confronted with studies of continu- 
ous Dyson-Schwinger equations (DSE). Both nonperturbative 
approaches have their own limitations. Whereas the lattice 
approach is affected by finite-size and discretization errors, 
the DSE approach requires truncations of an infinite tower of 
equations and those truncations are difficult to control. There- 
fore, a comparison of results is eventually able to provide 
more confidence about the consistence of both approaches. 

Starting with the work by L. von Smekal et a/. HH, DSE 
studies in recent years |Hl H, 0, S 111 [H have shown evi- 
dence for an intertwined infrared power behavior of the gluon 
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and ghost dressing functions 
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respectively, with the same 

value A. « 0.59 iH 113 ■ 

Thus, the 

gluon propagator D{q^) would be vanishing in the infrared 
in close connection with a diverging ghost propagator G{q^). 
This infrared behavior is closely related to gluon and quark 
confinement in accordance with the Gribov-Zwanziger hori- 
zon condition [13. , 141 and the Kugo-Ojima criterion [ is]]. 

As a by-product, a nonperturbative determination of the 
running coupling OLs{q^) in a momentum subtraction (MOM) 
scheme can be obtained. In fact, under the condition that the 
ghost-gluon-vertex renormalization function Zi(/j^) is finite 
and constant (see dill and d for a recent SU{2) lattice 
study) the corresponding running coupling is defined by 



(2) 



This together with relation ([TJm'ovides a non-trivial fixed 
point of ttj in the infrared limit [11], which was also proposed 
by D. V. Shirkov on the basis of a perturbative analytic ap- 
proach (see [19] and references therein). 

For quenched SU{2) extensive lattice investigations of the 
gluon and ghost propagators in Landau gauge can be found in 
ll20ll . For SU{3) lattice computations of the gluon propagators 
were reported already in |21l p. Later, other groups 



focused on the ghost propagator 11251 1261I27I1 . too. In our study 
we first paid special attention to the Gribov copy problem 
[i28] . In the latter context we have also investigated spectral 
properties of the Faddeev-Popov operator i29ll . With increas- 
ing physical volume we have found the low-lying eigenmode 
spectrum becoming steeper at the lowest non-zero eigenval- 
ues. This hopefully closes the gap from the trivial zero eigen- 
values and leads to a non-vanishing spectral density at zero 
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in the thermodynamic limit which is required for an infrared 
diverging ghost propagator. Moreover, in L30i1 we have re- 
ported on a first SU (3) lattice computation of the ghost-gluon 
vertex at zero gluon momentum that lends confirmation (see 
below) for an almost constant ghost-gluon-vertex renormal- 
ization constant. This was shown earlier in [18] for the case 
of SU (2). 

Also investigations for the gluon and ghost propagator in 
full QCD have been reported 0, [Hill. We have extended 
our investigations to this case, too, using configurations gener- 
ated with Nf = 2 dynamical clover-improved Wilson fermions 
lf34l . The lattice field configurations have been provided to us 
by the QCDSF collaboration | 35','3^. 

Finally, we only mention the analysis of some interrelated 
confinement criteria: a check of reflection positivity viola- 
tion by the gluon propagator [37J and the computation of the 
Kugo-Ojima confinement parameter Issll . Our own checks of 
these issues (see |]J]) have led to similar observations and will 
be published elsewhere ll39tl . 



II. LANDAU GAUGE, LATTICE GLUON AND GHOST 
PROPAGATORS 

The SU{3) lattice gauge field configurations U = {Ux.fi} 
have been generated by a standard (Hybrid) Monte Carlo algo- 
rithm and then put into the Landau gauge by iteratively maxi- 
mizing the gauge functional 



X,IJ J 



X 1.1=1 



''^t/v.A- =&t4,//^l+^ (3) 



with gx G SU{3). In general there are numerous local max- 
ima (Gribov copies), each satisfying the lattice Landau gauge 
condition 

{d,%){x)=Y,{%{x + ^/2)-^A,{x-ij/2))=0 (4) 



for the gauge transformed lattice potential 
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To explore to what extent this ambiguity has a significant in- 
fluence on gauge dependent observables, for lattice sizes up to 
24"* we have gauge fixed each thermalized configuration a cer- 
tain number of times (several tens depending on the lattice size 
and coupling constant) employing the over- relaxation algo- 
rithm and starting always from purely random gauge copies. 
Then for each configuration U, we have selected the first (f c) 
and the best (be) gauge copy (best with respect to the func- 
tional value) for subsequent measurements. For details we 
refer to lf28ll . On lattices sizes larger than 24^ we have re- 
stricted ourselves only to one copy per thermalized configura- 
tion (only f c copies). On those lattices we have applied also 
the Fourier accelerated method for fixing U to the Landau 
gauge. 



The momentum space gluon propagator is the cor- 

relator of two Fourier transforms A^(A:) of the lattice potential 
%^(x +/j/2) (the symbol is dropped from now on) 
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where q denotes the "physical" momentum 



qy,{k^,) 



( T^kn 



■ sm - 



(6) 



(7) 



related to the integer valued lattice momentum kjj G 
{~Lfj/2,Lfj/2] for the linear lattice extensions L^,/j — 1, ... ,4. 
According to Ref. ll22ll . a subset of admissible lattice mo- 
menta k has been chosen for the final analysis of the gluon 
(and ghost) propagator, although the Fast Fourier Transform 
algorithm easily provides us with all lattice momenta. 

The ghost propagator is derived from the Faddeev-Popov 
(F-P) operator, the Hessian with respect to of the gauge 
functional given in Eq. (O. It can be written in terms of the 
(gauge fixed) link variables Ux,fj as 



xy ~ l^^x.ii^x.y "x,iJ^x+iJ.y '^x./j^x-p.y 



(8) 



with A 
B 
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and r" . fl = 1 , . . . , 8 being the (hermitian) generators of the 
su(3) Lie algebra satisfying Tr [T"T''] = 5"'' /2. The ghost 
propagator is then determined by inverting the F-P operator 
M 



1 , 
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Following Refs. ll40ll4Tll we have used the conjugate gradient 
(CG) algorithm to invert M on a plane wave with color 
and position components "^"{x) = 5"' exp(27t/A;-x). In fact, 
we applied a pre-conditioned CG algorithm (PCG) to solve 
M«^(|)^(3;) = As the pre-conditioning matrix we used 

the inverse Laplacian operator ' with a diagonal color sub- 
structure. This has significantly reduced the required amount 
of computing time (for details see ll28ll ). 



III. PROPAGATORS: QUENCHED AND FULL QCD 
RESULTS 

For the pure SU (3) Yang-Mills case we have generated the 
gauge fields using the standard Wilson plaquette action for 
bare coupHng values P = 6/§q — {5.7,5.8,6.0,6.2}. Study- 
ing very large lattices (up to lattice size 56^) at the smallest 
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FIG. 1: The dressing functions for the gluon propagator Z{q^) = q^D{q^) vs. for quenched QCD (l.h.s.) and full QCD (r.h.s.), both 
measured on f c gauge copies. To illustrate the unquenching effect some quenched QCD data (P = 6.0, 32"* and 48"*) are shown in the right 
figure, too. 
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FIG. 2: As Fig.[T] here for the dressing function of the ghost propagator Jicp') = q^G{q^). The line (l.h.s.) represents a fit to the data for 
momenta lower than qf using the power law given in relation ([ij. 



10.0 



5.0 - 



0.0 - 




0.01 



0.1 



1 

g2 [GeV^ 



10 



100 



FIG. 3: The gluon propagator D{q^) vs. q^ for quenched QCD, 
measured on f c gauge copies. 



possible p value we tried to probe the infrared Umit, but with 
the reservation that the results can be affected by lattice dis- 
cretization artifacts. To match our lattice results to physical 
units we used the parametrization of the lattice spacing a de- 
pending on p as determined in [42]. For the full QCD case 
we have used Hybrid Monte Carlo generated QCDSF gauge 
field configurations produced with Nf — 2 dynamical flavors 
of clover-improved Wilson fermions. The relevant parame- 
ter sets 1.35. .36.1 are listed in Table |I] As lattice size we used 
always 24^ x 48. Note that the asymmetry demands some cau- 
tion in the very infrared. 

The gluon dressing function Z{q^) and the ghost one J{q^) 
are shown versus in physical units in Fig. [T] and Fig. |2] 
respectively, always for the f c copies only. The data represent 
the current status of our quenched and full QCD calculations. 
In particular, for the quenched case we have included very 
new results for the gluon propagator obtained at p = 5.7 using 
the lattice sizes 48^* and 56^*. All dressing functions have been 
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TABLE I: The parameter sets p, K etc. for the configurations used in 
our full QCD investigation (courtesy of the QCDSF collaboration). 



p 


K 




ma 


a[GeV"'] 


# conf 


5.29 


0.13550 


0.136410(09) 


0.0246 


2.197 


60 


5.29 


0.13590 


0.136410(09) 


0.0138 


2.324 


55 


5.25 


0.13575 


0.136250 (07) 


0.0135 


2.183 


60 



renormalized separately for each p such that they equal unity 
SLtq = 4 GeV. 

At a first glance our data for the gluon dressing function 
seem to be qualitatively in agreement with an infrared sup- 
pression, whereas the ghost dressing function looks compat- 
ible with an infrared singularity. The unquenching effect is 
clearly visible for the gluon propagator, whereas the ghost 
propagator is almost unaffected by including the fermionic 
feed-back in the functional measure of the gluons. This does 
not come unexpected since the ghost fields do not directly 
couple to the fermion fields. The non-perturbative peak of the 
gluon dressing function at ~ 1 GeV becomes softer as the 
quark mass is decreasing. It should be recalled that the gluon 
peak is wholly removed when center vortices are removed 
from the gluon field [43]. The corresponding effect of dy- 
namical quarks has been observed also in other lattice compu- 
tations of the gluon propagator using dynamical AsqTad im- 
proved staggered quarks 1 31] and is expected from studies of 
the unquenched ghost and gluon propagators within the DSE 
approach I 44i i45i [46(1 . We refer also to lattice studies with 
dynamical Kogut-Susskind and Wilson fermions reported in 
ll32i[33ll and to the contribution given by S. Furui during this 
meeting llssll . 

For the gluon as well as for the ghost dressing functions 
in the quenched case, with somewhat improved data in com- 
parison to [28,1 . we have tried fits of a power- like infrared be- 
havior (see iM)- For the ghost dressing function the result is 
indicated in the l.h.s. of Fig. |2] The resulting exponent turned 
out to be quite stable against variations of the upper end of the 
fit interval qj. We have found X = 0.20(1), i.e. much smaller 
than expected from the DSE approach. A similar observa- 
tion was made for the gluon exponent. However, in order to 
check whether the gluon propagator really vanishes in the in- 
frared — in accordance with the Gribov-Zwanziger horizon 
condition — we show the gluon propagator itself in Fig. |3] 
With our new quenched SU (3) data obtained at p = 5.7 on a 
56''^ lattice we now find a first indication for having reached 
a maximum, opening the opportunity for a decreasing gluon 
propagator for even lower momenta. Our results do not seem 
to expose dramatic finite-size effects, but obviously one needs 
much larger lattices in order to reach momenta well below 
100 MeV. Therefore, it is premature to attempt fits of the in- 
frared exponents with the hope to reproduce the DSE results 
(see Eq. ([T])). This certainly also holds for the ghost propaga- 
tor, i.e. the X value quoted above cannot be taken seriously. 



IV. GHOST-GLUON- VERTEX FUNCTION AND THE 
RUNNING COUPLING 

In Fig.|4]we show the running coupling according to Eq. (|2]). 
Surprisingly, the coupling is seen to decrease below ~ 
0.3 GeV^ rather than to approach the predicted non-trivial in- 
frared fix-point monotonously from below. The same hap- 
pens as well for the quenched as for the unquenched case, ir- 
respective of the clearly visible unquenching effects. Whether 
strong finite-size effects in the ghost propagator may change 
the present tendency remains to be seen in future. From our 
checks of finite-size effects (see below) we cannot derive ar- 
guments that this will be the case. It is worth mentioning 
here that lattice computations of the running coupling from 
other vertex functions have provided quite similar results (see 
ll47i |48[| for the three-gluon vertex and |49] for the quark- 
gluon vertex). 

The definition of the running coupling relies on the assump- 
tion of a constant ghost-gluon-vertexrenormalization function 
Zi {q^). A recent lattice investigation of this function defined 
at vanishing gluon momentum for the SU{2) case (l^ sup- 
ports that Zi (q^) ^ I at least for momenta larger than 1 GeV. 
We have performed an analogous study for Zi {q^) in the case 
of SU (3) gluodynamics and for full QCD. Our results are pre- 
sented in Fig.|5] There is a slight variation visible in the inter- 
val 0.3 GeV^ < < 5 GeV^. However, this weak deviation 
from being constant would not have a dramatic influence on 
the running coupling. 

The left-most two data points falling below unity on the 
rh.s. of Fig.|5]correspond to the lowest on-axis momentum on 
the (asymmetric) 24^ x 48 lattice. We believe that this devi- 
ation is due to the asymmetry of the lattice, because we have 
seen a similar effect (not shown) in our quenched data, too. 
Simulations on larger (symmetric) lattices will enable us to 
sharpen our conclusions within the near future. 



V. SYSTEMATIC LATTICE AND GRIBOV COPY 
EFFECTS 

It is a common problem of all lattice computations that one 
has to check the sensitivity of the results with respect to finite- 
size effects, to the boundary conditions adopted, to an even- 
tual lattice asymmetry, and to the given lattice discretization. 
We discuss three of these checks in the following. Note that 
we have employed some cuts on the list of momenta used in 
order to minimize effects of the lattice discretization from the 
beginning (see ll22ll.lfTll). 

In Fig.|6]we demonstrate how the finite lattice size can influ- 
ence the propagator results. For several fixed p values, i.e. for 
fixed lattice spacings we check the dependence on the varying 
volume. For smallest lattices (at smallest lattice spacings for 
P — 6.2) we find very strong finite-size effects in both cases as 
one has expected. But our results show that at larger volumes 
(at lower p values) the lattice size dependence becomes less 
dramatic. 

As a part of the project we also performed simulations at 
P = 6.0 using various asymmetric lattices. In Fig. [7] we show 
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FIG. 4: The momentum dependence of the running couphng as{q^) for quenched QCD (l.h.s.) and full QCD (r.h.s.)i, measured on f c gauge 
copies. For comparison, on the right hand side selected quenched QCD data for (3 = 6.0, 32^ are shown. 1- and 2-loop fits to CXj(^^) are drawn 
with dashed and solid lines, respectively. 
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FIG. 5: The inverse ghost-gluon-vertex renormalization function Zj ' (q^), measured on f c gauge copies for quenched QCD (l.h.s.) and full 
QCD (r.h.s.). 



gluon (l.h.s.) and ghost (r.h.s.) dressing function data obtained 
for lattice sizes 16^ x 128 and 24^ x 128 in comparison with 
the symmetric lattice 48"*^. We see large systematic effects at 
low momenta due to the asymmetry, in particular, for the low- 
est on-axis momenta along the elongated time direction. For a 
more detailed study of this effect see 1 1]. Similar observations 
have been made from 3d SU (2) investigations in |50]. At this 
workshop P. Silva has reported on the possibility to use the 
asymmetry effect for an extrapolation towards the large vol- 
ume limit lIsTl I52I1 . It remains to be seen, whether such an 
extrapolation really gives stable results. 

Next, in Fig. [8] we present a check of the dependence on 
the lattice discretization for the ghost and the gluon dress- 
ing functions in quenched QCD. Whereas the lattice volumes 
are approximately the same, the lattice spacing changes from 
a ~ 0.136 fm (p = 5.8) to a ~ 0.093 fm (p = 6.0). The 
comparison shows that lattice discretization effects are almost 
negligible (for the ghost) or quite moderate (for the gluon) in 
the given momentum interval for the set of (preselected) lat- 
tice momenta. 



Finally, in Fig.|9]we illustrate the effect of the Gribov copies 
for periodic gauge transformations. We have plotted some f c 
- to - be ratios of the ghost and gluon dressing functions. 
For the gluon propagator there is no influence revealing itself 
on top of the statistical noise. On the contrary, for the ghost 
propagator the Gribov problem can cause (9(5%) deviations 
in the low momentum region (^ < 1 GeV). For better gauge 
copies the ghost dressing function becomes less singular in 
the infrared. A closer inspection of the data for the ghost 
propagator indicates that the influence of Gribov copies be- 
comes weaker for increasing physical volume. This probably 
corroborates a recent claim by Zwanziger telling that in the 
infinite volume limit averaging over gauge copies in the Gri- 
bov region should lead to the same result as averaging over 
copies restricted to the fundamental modular region 15 311 . In 
[54] we have investigated for SU{2) gluodynamics the effect 
of admitting a non-periodic extension (modulo Z(2) flips) of 
the usually periodic gauge transformations. Also in this case 
we found that the effect on the Greens functions of the bias to- 
wards a particular gauge copy fades away with increasing lat- 
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FIG. 6: Demonstration of finite-size effects of the dressing functions for tfie ghost (upper row) and for the gluon (lower row) in quenched 
QCD. 
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FIG. 7: Lattice asymmetry effects for the gluon (l.h.s.) and ghost (r.h.s.) propagators, respectively, in quenched QCD (p = 6.0, only f c gauge 
copies). Note that we have dropped data points at even lower momentum (on-axis in the elongated time direction). Our intention here is to 
concentrate on the momentum range where we can compare with data obtained on symmetric lattices. 



tice volume. Anyway, we see that Gribov copy effects have to 
be studied properly in the infrared limit before one can come 
to final conclusions. 



VI. SUMMARY 

We have studied the low momentum region of QCD in the 
Landau gauge using Monte Carlo simulations with the Wilson 
plaquette action. For the quenched case we could afford to 



simulate on lattice sizes ranging from 8'* up to 56^ using bare 
couplings constants p in the interval [5.7, ... ,6.2]. In this way 
we have presently reached momenta down to ^ ~ 100 MeV. 
In order to assess the effect of virtual quarks we have carried 
out an analogous investigation for full QCD with two flavors 
of clover-improved dynamical Wilson fermions for three dif- 
ferent quark masses. For the time being we have evaluated 
only lattices of size 24^ x 48. Our data presented here refer 
to an arbitrary (the "first") gauge copy as obtained from over- 
relaxation or Fourier accelerated gauge fixing. 
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FIG. 8: Lattice discretization effects of the ghost (l.h.s.) and gluon (r.h.s.) dressing functions, respectively, in the quenched case. Both dressing 
functions are shown for two cases (P = 5.8, 16^ and (3 = 6.0,24'*) corresponding to approximately the same lattice volume ~ (2.2 fm)^. 
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FIG. 9: The ratios Z'^'^/Z'"^ for the gluon dressing functions (l.h.s.) and y'^^/i'"^ for the ghost dressing functions (r.h.s) determined on first (f c) 
and best (be) gauge copies, respectively, as function of the momentum . 



Towards the infrared momentum region, the gluon dressing 
functions in quenched as well as in full QCD were shown to 
decrease, while the ghost dressing functions turned out to rise. 
However, the interrelated power laws predicted by the infinite- 
volume DSE approach could not (yet) be confirmed on the 
basis of our data. Our new data in the quenched case obtained 
for lattice sizes 48'* and 56^ at p = 5.7 demonstrate that the 
gluon propagator flattens for < 0.1 GeV^ leaving open the 
possibility for its decrease at even lower momenta. 

From the present data, the running coupling (Xs{t') in the 
momentum subtraction scheme (based on the ghost-gluon ver- 
tex) does not seem to approach the expected finite infrared 
fixed point monotonously. It was rather seen to decrease for 
lower momenta after passing a turnover at ~ 0.3 GeV^ for 
quenched {Nf = 0) as well as for full QCD {Nj = 2). 

Unquenching effects have been clearly identified for the 
gluon propagator, whereas the ghost propagator was almost 
unchanged. This is in one-to-one correspondence with what 
has been found in the Dyson-Schwinger equation approach. 
However, the puzzle of the existence of a non-trivial infrared 
fixed point in the infinite volume limit remains unsolved. 



We have studied lattice effects as far as the dependence on 
the finite lattice size, on the lattice spacing and on the lattice 
asymmetry are concerned. Within the parameter range under 
study the first two problems seem to be under control, whereas 
infinite volume extrapolations based on the strong dependence 
on the lattice asymmetry might be worth to be further studied. 

Concerning the effect of Gribov copies we have seen a quite 
strong influence in the infrared region on the ghost propaga- 
tor which becomes less singular when better gauge copies are 
taken. We have found some indications that the Gribov effect 
weakens as the volume increases. 
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